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ABSTRACT 

Context. Asymptotic giant branch (AGB) stars are in one of the latest evolutionary stages of low to intermediate-mass stars. Their 
vigorous mass loss has a significant effect on the stellar evolution, and is a significant source of heavy elements and dust grains for 
the interstellar medium. The mass-loss rate can be well traced by carbon monoxide (CO) line emission. 

Aims. We present new Herschel HIFI and IRAM 30m telescope CO line data for a sample of 53 galactic AGB stars. The lines cover 
a fairly large range of excitation energy from the 7 = 1 —> 0 line to the 7 = 9 —> 8 line, and even the 7 = 14 —> 13 line in a few cases. 
We perform radiative transfer modelling for 38 of these sources to estimate their mass-loss rates. 

Methods. We used a radiative transfer code based on the Monte Carlo method to model the CO line emission. We assume spherically 
symmetric circumstellar envelopes that are formed by a constant mass-loss rate through a smoothly accelerating wind. 

Results. We find models that are consistent across a broad range of CO lines for most of the stars in our sample, i.e., a large number 
of the circumstellar envelopes can be described with a constant mass-loss rate. We also find that an accelerating wind is required to 
fit, in particular, the higher-J lines and that a velocity law will have a significant effect on the model line intensities. The results cover 
a wide range of mass-loss rates (~ 10“* to 2 x 10“^ Mq yr“') and gas expansion velocities (2 to 21.5 km s^*), and include M-, S-, and 
C-type AGB stars. Our results generally agree with those of earlier studies, although we tend to find slightly lower mass-loss rates by 
about 40%, on average. We also present “bonus” lines detected during our CO observations. 

Key words. Stars: AGB and post-AGB - circumstellar matter - stars: mass loss - stars: evolution 


1. Introduction 

Towards the end of their lives, low and intermediate mass stars 
(with masses ~ 0.8-8 Mq) will exhaust their supply of He and 
cease fusion reactions in their cores, leaving a quiescent C/O 
core with H and He fusion reactions only taking place in thin 
shells surrounding the core. This evolutionary phase is known as 
the asymptotic giant branch (AGB) ( |Herwig|2005| l. 

AGB stars are also a significant source of heavy elements in 
the universe. It is thought that about half of all elements heav¬ 
ier than Ee originate in AGB stars through the i-process of slow 
neutron capture ( |Herwig|2005| l. It is during the AGB phase that 
this enriched material is brought to the surface. At the same time, 
the star experiences vigorous mass loss, ejecting matter to form 
a circumstellar envelope (CSE) around the star. Molecules and 
dust grains form in the CSE, and will eventually chemically en¬ 
rich the interstellar medium (ISM). 


* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im¬ 
portant participation from NASA. 

** Based on observations carried out with the IRAM 30m Telescope. 
IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 


It is believed that AGB stars begin their lives on the AGB as 
oxygen-rich M-type stars, and eventually some of these, those 
with masses in the range ~ 1.5^ Mq ( Herwig|2005| l, will tran¬ 
sition into carbon-rich C stars. With a C/O ratio close to 1, S 
stars are believed to occupy the evolutionary phase between M 
and C stars. The lowest mass AGB stars (< 1 Mq) do not trans¬ 
form into C stars because they do not undergo a third dredge-up 
event. The highest mass AGB stars (> 4 Mq) also do not end 
their lives as C stars due to hot bottom burning (HBB), unless 
the mass loss quenches the HBB process, leaving time for the 
star to evolve into a C star before leaving the AGB. 

Radiative transfer modelling of circumstellar CO radio 
lines has long been used to derive the mass-loss rates of 
AGB stars (|Morris||1987[ |Kastner||1992[ [Justtanont et al.||1994t 


Groenewegen|| 1998 Schoier & Oiofsson||2001[ Olofsson et al. 


Decin et al.| |2006t [Ramstedt et al.| |2b09[ |De Be^ 


2002 


et al. 


|2010l. These data were almost exclusively obtained with 


ground-based telescopes. The Herschel Space Observatory al¬ 
lowed observations of higher energy lines than possible from 
ground-based telescopes. This has led to studies that model 
molecular emission (not only that of CO) in more detail and over 
a wide range of energies as in Schoier et al. (2011 1 , Khouri et al. 
(|2014|l, and|Danilovich et al.|(2014[). However, each of those pa- 
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pers deals with only one or a handful of stars. The new data have 
not yet been applied to mass-loss rate determinations for larger 
samples of stars. 

In this paper we present new data obtained as part of the 
SUCCESS project (SUbmilimetre Catalogue of Circumstellar 
Envelopes with HerSchellY\W\, Teyssier et al.|20lT i, a CO sur¬ 
vey of a large sample of AGB and post-AGB stars across all 
three chemical types. HerschellWWl was used to observe the CO 
7 = 5 —> 4 and 7 = 9 —> 8 lines in 53 AGB stars. Of these, six 
stars were also observed in the CO J - 14 —» 13 line and 29 
stars were also observed in the CO 7 = 1 —» 0 and 7 = 2 —> 1 
lines using the IRAM 30m telescope. Of these stars we present 
radiative transfer models for 38 objects, supplementing our new 
high-resolution observations with archival data where possible. 
Inclusion of the high-7 lines from HerschellHWl allows us to 
better constrain our models over a broad range of temperatures, 
and hence produce models that better represent the CSE over a 
large radial range. 


2. Sample and observations 

SUCCESS is a HerschellHWl Guaranteed Time project 


(Teyssier et al. 2011 1 which observed 74 AGB and post-AGB 
stars in the CO 7 = 5 —> 4 and 7 = 9^8 lines, and 10 of 
those objects also in the CO 7 = 14 —» 13 line. Excluding post- 
AGB objects and extreme OH/IR stars (the latter presented in 
lJusttanont et al.|[2013[ ), the sample comprised of 53 AGB stars 
observed in the former two lines and six observed in all three 
lines. The sample was primarily built based on the evolved star 
lists considered in the st udies by [Knapp et al.| ([T998| l, |Schoier| 
& Olofsson (|2001|l, and Castro-Carrizo et al. (|2010|l, comple 


mented by some of the brightest line calibrators regularly used at 
the APEX telescope ( |Gusten et al.|2008l l. The selection criterion 
was related to the intensity of the CO 7 = 2 —> 1 as observed at 
the IRAM 30m telescope, and aimed at sources exhibiting peak 
intensities in excess of 1.5 K in this line. 

A reduced sample of 29 objects was also observed with the 
IRAM 30m telescope in the CO 7 = 1 —> 0 and 7 = 2^1 lines. 
The objective was to obtain well-calibrated data observed with 
state-of-the-art receivers, especially for those sources where the 
literature spectra dated from several decades ago and/or had been 
obtained from a varied set of facilities. 

The full sample is summarised in Table including the pul¬ 
sation periods and the systemic velocities. A summary of all our 
new observations and integrated line intensities is given in Table 
1^ We have indicated with an * those lines which have some ISM 
contamination that was corrected for when calculating the inte¬ 
grated intensity. In these cases, the contamination tended to oc¬ 
cur only on one side of the line and we calculated the intensity by 
integrating from the centre to the edge of the non-contaminated 
side and then doubling this result. The beam-widths for each line 
and telescope are given in Table 

The new data for stars we have modelled are shown in 
Eigures |A.1| |A.2| and |A.3| The new observations for stars we 
have not modelled are shown in Figures [BTT]|B.2| and B.3 


2.1. IRAM observations 

The complementary observations at the IRAM 30m telescope 
were obtained in December 2013 respectively using the E090 
and E230 EMIR receivers and the FTS backend to cover the 
CO 7 = 1 —» 0 and 7 = 2 —> 1 lines. We used the Wobbler 
Switching mode with a throw of 4'. Owing to the large instanta¬ 
neous bandwidth offered by the EMIR receivers, several bonus 


lines from '^CO, CN, SiS, and HC 3 N were simultaneously ob¬ 
served and are presented in Section 2.4 Data were processed 
with the GILDAS/Clas^ software and converted into the Tinb 
scale assuming main beam efhciencies of 0.78 and 0.58 at 3 and 
1 mm respectively. We have corrected for the antenna elevation 
gain, which accounts for up to 25% line intensity loss in our data. 
Finally, we re-adjusted our line intensities based on the monitor¬ 
ing of reference spectra of the CO lines in IRC-n 10216 (with 
hducial line peak intensities of 24.3 K and 54.5 K for the CO 
7 = 1 —> 0 and 7 = 2 —> 1 lines, respectively), resulting in cor¬ 
rections between 0 and 20 % depending on the line and the ob¬ 
serving day. The achieved noise rms are in the range 30-150 mK 
and 55-360 mK at 3 and 1 mm, respectively, and per native ve¬ 
locity resolution element (0.5 and 0.25 km s“\ respectively). 


2.2. HIFI observations 


The SUCCESS sample was observed between July 2010 and 
April 2012 using the HIFI instrument (de Graauw et al. 2010| l 
aboard the Herschel Space Observatory (Pilbratt et al. 2010 1 . 
The data were obtained with the double beam switching mode, 
using reference positions separated by 3' from the target posi¬ 
tion ( Roelfsema et al.|[2012 l. The spectra were sampled on the 
Wide Band Acousto-Optical Spectrometer (WBS), offering a na¬ 
tive resolution of 1.1 MHz (0.6 kms * at the lowest observed 
frequency). The selected frequency tunings were optimised for 
the targeted CO lines, but bonus lines of CN and SiO were also 
obtained in the instantaneous bandwidth of 4 GHz (2.4 GHz for 
the CO 7 = 14 —> 13 line) provided by the WBS (see Sect. 
|2.4| l. The full list of Herschel observation identihers (ObsIDs) is 
given in Table |B.l| 

The HIFI data have been processed using HIPE l^and cal¬ 
ibrated in the scale. On top of that, the CO 7 = 14 —> 13 
data have been corrected for the so-called electrical standing 
waves using the doHebCorrection task (Kester et al. 2014| l. A 
sideband gain ratio different from unity has been used for the 
CO 7 = 5^4 observations ( [Higgins et al.||2014) l. Finally, all 
our data have been converted to the using the revised main 
beam efficiencies (Mueller et al., 201^. Likewise, we refer to 
this technical note for the beam size assumed in our modelling 
(Table [^. The achieved 1-cr noise rms was 15 mK (T^b) for a 
smoothed resolution channel of 3 km s * for both CO 7 = 5^4 
and 7 = 9 —> 8 , and 44 mK for CO 7 = 14 —> 13. 


2.3. Supplementary data 

To better constrain our models, we used previously observed 
low- and intermediate-7 CO lines from a variety of sources. A 
summary of the supplementary observations is given in Table 

01 

The telescopes and their corresponding beam widths for dif¬ 
ferent frequencies are listed Table covering all the new and 
archival observations used in this paper. 

As well as archival line data, we have used results from the 
APEX Pointing Catalogue, which can be found onlin^ 


‘ http://www.iram.fr/IRAMFR/GILDAS 

^ HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS, and SPIRE consortia. 

^ http://herschel.esac.esa.int/twiki/pub/Public/HifiCalibrationWeb/ 
HifiBeamReleaseNote_Sep2014.pdf 

http://www.apex-telescope.org/observing/pointing/spectra/ 
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Table 1. Basic stellar parameters for the AGB stars in the SUCCESS sample 


Star 

Alternate name 

RA 

(J2000) 

Dec 

(J2000) 

Variability 

Period 

[days] 

Ulsr 

[km s“*] 

C stars 

RScl 

IRC-30015 

01 26 58.09 

-32 32 35.5 

SRb 

370 

-19 

V701 Cas 

AFGL 482 

03 23 36.57 

+70 27 07.5 

M 

567 

-14 

V384 Per 

IRC-H50096 

03 26 29.53 

+47 31 50.2 

M 

535 

-16 

GY Cam 

IRC+60144 

04 35 17.54 

+62 16 23.8 

SR 


-47 

R Lep 

IRC-10080 

04 59 36.35 

-14 48 22.5 

M 

445 

11 

V1259 Ori 

AFGL 865 

06 03 59.84 

+07 25 54.4 

M 

696 

42 

IRAS 06192-H4657 


06 22 58.52 

+46 55 34.9 

7 


-23 

UU Aur 

IRC +40158 

06 36 32.84 

+38 26 43.8 

SRb 

235 

8 

V688 Mon 

AFGL 971 

06 36 54.24 

+03 25 28.7 

M 

653 

3 

AI Vol 

IRAS 07454-7112 

07 45 02.80 

-71 19 43.2 

M 

511 

-39 

UHya 

IRC-10242 

10 37 33.27 

-13 23 04.4 

SRb 

183 

-32 

XTrA 

IRAS 15094-6953 

15 14 19.18 

-70 04 46.1 

SR 

361 

-2 

II Lup 

IRAS 15194-5115 

15 23 04.91 

-51 25 59.0 

M 

575 

-15.5 

VCrB 

IRC+40273 

15 49 31.31 

+39 34 17.9 

M 

358 

-99 

V82I Her 

IRC+20370 

18 41 54.39 

+ 17 41 08.5 

M 

524 

-0.5 

V Aql 

IRC-10486 

19 04 24.16 

-05 41 05.4 

SRb 

407 

54 

VI968 Cyg 

AFGL 2494 

20 01 08.51 

+40 55 40.2 

M 

783 

28 

RV Aqr 

IRC+00499 

21 05 51.68 

-00 12 40.3 

M 

453 

1 

TXPsc 

IRC+00532 

23 46 23.52 

+03 29 12.5 

LB 


13 

S stars 

R And 

IRC+40009 

00 24 01.95 

+38 34 37.3 

M 

409 

-16 

SCas 

IRC+70024 

01 19 41.97 

+72 36 39.3 

M 

612 

-30 

W And 

IRC+40037 

02 17 32.96 

+44 18 17.8 

M 

397 

-35 

R Gem 

IRC+2017I 

07 07 21.27 

+22 42 12.7 

M 

370 

-60 

YLyn 

IRC+50180 

07 28 11.61 

+45 59 26.2 

SRc 

110 

-0.5 

RS Cnc 

IRC+30209 

09 10 38.80 

+30 57 47.3 

SRb 

242 

6.5 

RCyg 

IRC+5030I 

19 36 49.38 

+50 11 59.5 

M 

426 

-17 

n Gru 

AFGL 4289 

22 22 44.21 

-45 56 52.6 

SRb 

196 

-12 

M stars 

KU And 

IRC+40004 

00 06 52.94 

+43 05 00.0 

M 

720 

-22 

V370 And 

IRC+50049 

01 58 44.33 

+45 26 06.9 

SRb 

228 

-2 

AFGL 292 

IRC+I0025 

02 02 38.63 

+07 40 36.5 

7 


23 

RHor 

IRAS 02522 -5005 

02 53 52.77 

-49 53 22.7 

M 

408 

37 

NV Aur 

IRC+50137 

05 11 19.44 

+52 52 33.2 

M 

635 

2 

BX Cam 

IRC+70066 

05 46 44.10 

+69 58 25.2 

M 


-2 

GXMon 

IRC+I0I43 

06 52 46.91 

+08 25 19.0 

M 

527 

-9 

L 2 Pup 

IRAS 07120 -4433 

07 13 32.32 

-44 38 23.1 

SRb 

141 

33.5 

SCMi 

IRC+I0167 

07 32 43.07 

+08 19 05.2 

M 

333 

51 

RLMi 

IRC+302I5 

09 45 34.28 

+34 30 42.8 

M 

372 

0 

R Leo 

IRC+I02I5 

09 47 33.49 

+ 11 25 43.7 

M 

310 

0 

RCrt 

IRC-20222 

11 00 33.85 

-18 19 29.6 

SRb 

160 

11 

BK Vir 

IRC+00220 

12 30 21.01 

+04 24 59.2 

SRb 

140 

17 

YUMa 

IRC+60220 

12 40 21.28 

+55 50 47.6 

SRb 

168 

19 

RT Vir 

IRC+I0262 

13 02 37.98 

+05 11 08.4 

SRb 

158 

18 

SW Vir 

IRC+00230 

13 14 04.39 

-02 48 25.2 

SRb 

146 

-11 

RHya 

IRC-20254 

13 29 42.78 

-23 16 52.8 

M 

380 

-10 

RXBoo 

IRC+30257 

14 24 11.63 

+25 42 13.4 

SRb 

158 

2 

SCrB 

IRC+30272 

15 21 23.96 

+31 22 02.6 

M 

360 

0 

XHer 

IRC+50248 

16 02 39.17 

+47 14 25.3 

SRb 

102 

-73 

VIIII Oph 

IRC+I0365 

18 37 19.26 

+ 10 25 42.2 

M 


-32 

RR Aql 

IRC+00458 

19 57 36.06 

-01 53 11.3 

M 

395 

28 

VI943 Sgr 

IRC-30425 

20 06 55.24 

-27 13 29.8 

SRb 

330 

-15 

VI300 Aql 

IRC-I0529 

20 10 27.87 

-06 16 13.6 

M 

680 

-18 

TCep 

IRC+70168 

21 09 31.78 

+68 29 27.2 

M 

388 

-2 

EP Aqr 

IRC+00509 

21 46 31.85 

-02 12 45.9 

SRb 

55 

-34 


Notes. An ellipsis (...) indicates an unknown property. Variability and period information was obtained from the International Variable Star Index 
(VSX) database. 


For our SED models we primarily used photometry from 2.4. Bonus lines 
IRAS and 2MASS ( |H0g et al.|2000l l. See Sect.j^for more de¬ 
tails. In the course of our new HIFI and IRAM observations, we also 

acquired some “bonus” line spectra for molecules that were ob¬ 
servable within our target frequency ranges. Our HIFI observa- 
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Table 3. Telescope parameters for all lines referred to in this 
paper 


Transition 

Frequency 

[GHz] 

Telescope 

e 

["] 

CO (1 ^ 0) 

115.271 

IRAM 

21.4 



NRAO 

55 



OSO 

33 



SEST 

45 

CO (2 ^ 1) 

230.538 

APEX 

27 



CSO 

30 



IRAM 

10.7 



JCMT 

21 



SEST 

23 

CO (3 ^ 2) 

345.796 

APEX 

18 



CSO 

20 



JCMT 

14 



SEST 

15 

CO (4 ^ 3) 

461.041 

APEX 

14 



CSO 

15.5 



JCMT 

12 

CO (5 ^ 4) 

576.268 

HIFI 

36.1 

CO (6 ^ 5) 

691.473 

CSO 

10.3 



HIFI 

30.4 



JCMT 

8 

CO (7 ^ 6) 

806.652 

APEX 

7.7 

CO (9 ^ 8) 

1036.912 

HIFI 

20.1 

CO (10^ 9) 

1151.985 

HIFI 

18.2 

CO (14^ 13) 

1611.794 

HIFI 

12.9 

CO (16^ 15) 

1841.345 

HIFI 

11.5 


Notes. APEX is the Atacama Pathfinder Experiment; CSO is the 
Caltech Submillimeter Observatory; IRAM refers to the 30m telescope 
at the Institut de Radioastronomie Millimetrique; JCMT is the James 
Clerk Maxwell Telescope; HIFI is the Heterodyne Instrument for the 
Far-Infrared aboard Herschei, NRAO refers to the 12 m telescope at the 
(US) National Radio Astronomy Observatory; OSO is the 20 m tele¬ 
scope at the Onsala Space Observatory; SEST is the Swedish-ESO sub¬ 
millimetre telescope. 


tions covered the CN ( 59/2 ^ 47 / 2 ) and (5n/2 ^ 49 / 2 ) line 
groups, which were detected in a handful of C stars, and SiO 
(13 —> 12), which was detected mostly in M stars. Our IRAM 
observations covered the CN N - I ^ 0 and CN N - 2 ^ I 
line groups, which were detected in most C stars and one S star 
(S Cas). Also covered by IRAM were the '^CO (1 —» 0) line, 
which was detected in most observed sources, the SiS (6 ^ 5) 
line, which was detected in higher mass-loss rate sources, and 
the HC 3 N (12 ^ 11) line, which was detected in the higher 
mass-loss rate C stars. 

We will not be modelling these additional species but we 
discuss the bonus detections in greater detail in Sect, 
observations are plotted in Figures |B.5| |B.4| |B. 6 | and 
integrated intensities for the detected lines are listed in Table 


B.2 


B.7 


The 

The 


3. Radiative transfer modeiiing 

3.1. SED modelling 

We begin our radiative transfer modelling by estimating some 
key dust properties of each star in the sample. The SED mod¬ 
elling is performed using DUSTTp] a publicly available radia¬ 
tive transfer code ( |Ivezic & Elitzur||l997[ ). We found the best 
fit for each star using primarily 2MASS and IRAS photometric 

^ http://www.pa.uky.edu/~moshe/dusty/ 


observations. The distances taken from the literature and the lu¬ 
minosities were calculated from the period-luminosity relation 
of IGlass & E\^ ( fTMT] l: 

Mto/= 0.76-2.09 log P (1) 


The resulting effective temperature of the central black body, 
dust optical depth (given at 10 fim) and the inner radius of our 
model, based on the dust condensation temperature, are listed in 
Table g] 


As the DUSTY code has been widely used, some of the stars 
in our sample have already been modelled using the same meth¬ 
ods and in these cases we simply use the earlier results from 
Schoier et al.|(2007] l, [Schoier et al.| ( |2013| l, and [Ramstedt &| 


Olofssonl 20141. These are indicated with a t in Table |4| 


3.2. CO line modelling 
3.2.1. The circumstellar model 


In our radiative transfer modelling of the observed CO lines, 
we assumed a spherically symmetric CSE, formed through con¬ 
stant and isotropic mass loss with a smoothly accelerating wind. 
Some of our modelled stars, such as GY Cam and AEGL 292, 
show line profiles which we are not able to reproduce within our 
adopted circumstellar model. See Sect. 4.4 for further discussion 
of this. 

The radial gas velocity law used in our modelling is given by 


v{r) = Wmin + (Uoc 



( 2 ) 


where Umm - 3 kms * is the approximate sound speed at the 
dust condensation radius, Uoc is the observed terminal expansion 
velocity, is the dust condensation radius as calculated in the 
corresponding SED model, and y 6 is a parameter used to adjust 
the acceleration in the inner part of the envelope. In general, we 
assume j8 — 1.0 but adjust it for the stars where this gives a no¬ 
ticeable improvement in the fits to the line shapes. The j3 values 
used for each star are given in Table Eor some stars we use a 
constant expansion velocity, either because the line width is es¬ 
pecially low (such as for S CMi and L 2 Pup with Voc -2 km s ') 
or because the constant velocity model was a much better fit to 
the observations than a model with a velocity profile. The range 
of results is discussed in more detail in Sect. 14.41 We also em¬ 
ploy a constant turbulent gas velocity of 0.5 kms * for all the 
modelled stars. 

A similar function to Eq.|^is used for the dust and drift ve¬ 
locities, where the drift velocity is the difference between the 
dust and gas velocities. The terminal drift velocity is assumed to 
be (for details see |Kwok|1975| l: 


t^drifcco = (3) 

V Me 


where L* is the stellar luminosity, M is the mass-loss rate, Q is 
an efficiency factor assumed to be 0.03 ( |Ramstedt et al.|2008| l, 
and c is the speed of light. The drift velocity profile is calculated 
from Eq. I^with this terminal velocity and with - 1- For 
each star, the same and (3 values are used for the drift (and 
hence dust) velocity as for the corresponding gas velocity. Where 
the gas velocity has /Jgas = 0 , we assume that the drift velocity 
has jSdrift = 1- Using this formulation, we find drift velocities 
in the range 2 kms“^- 18 kms“* (with a median of 6 kms^'), 
which is below the sputtering-dominated threshold of 20 kms"* 
calculated by |Kwok| ( |1975| l. 
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3.2.2. The modelling approach 

We modelled the observed CO lines using a Monte Carlo pro¬ 


gram (MCP) which has been previously described in |Bernes| 


( 1979|l, 

Schoier & Olofsson ( 

2001 Jl, Schoier et al.| ^2002), 

Rams tec 

t et al. (2008 1 , and Danilovich et al. (2014 1 . The MCP 


code takes basic stellar and molecular parameters, and the re¬ 
sults of our SED model as input and calculates the molecular 
excitation by solving the statistical equilibrium equations using 
the Monte Carlo method. MCP also solves the energy balance 
equation to calculate the gas temperature as a function of radius 
throughout the CSE, 


rfTkin 

dr 


(2-2r) 



r dv\ 
2v{r) dr j 


Tkin(r) 


y- 1 

nn2kBv(r) 


(H-cm 


where Tkin is the kinetic temperature of the gas, is the hy¬ 
drogen number density, kg is Boltzmann’s constant, y is the adi¬ 
abatic index with y = f for Tkin < 350 K and y - ^ otherwise, 
and H and C are the sums of the heating and cooling terms, re¬ 
spectively. The adiabatic cooling is given by the hrst term on 
the right-hand side of the equation. Otherwise, the cooling terms 
include H 2 vibrational line cooling and CO rotational line cool¬ 
ing, and the heating terms include heating due to dust-gas colli¬ 
sions and photoelectric heating. CO line cooling can also act as 
heating in some circumstances (for more details, see |Schoier &| 
|01ofsson|2001| l. 

The most important heating term in the energy balance equa¬ 
tion, the dust-gas collision term, includes a number of assump¬ 
tions on the dust properties. These come in a multiplicative form 
and we combine them in the so called h-parameter dehned as 


/ T* \ / 2.0 gcm-3W 0.05 jum 

\0.0lj\ pd ag 


(5) 


where 'P is the dust-to-gas ratio, pd is the dust grain density and 
Qg is the dust grain radius, assuming spherical grains. We assume 
Pd - 2.2 g cm“^ for carbon dust and pd - 3.3 g cm~^ for silicate 
dust, and Og - 0.05 pm. The h parameter is one of the two free 
parameters adjusted (in practice we vary T') in our modelling to 
ht the observational data, the other being the mass-loss rate. The 
impact of the h-parameter is on the kinetic temperature structure 
of the circumstellar envelope, which tends to be close to a power 
law for part of the envelope. The molecular excitation analysis is 
complex and therefore the effect of changing h is not always eas¬ 
ily predictable, but the tuning of the h parameter hnally results in 
a kinetic temperature distribution that best matches the observed 
data, in particular the relative intensities of the different J lines. 
However, since we may be missing terms in the energy balance 
equation, or the included terms do not fully capture the physics 
involved, it is clear that the resulting T* of the best-ht model is 
not necessarily a good estimate of the circumstellar dust-to-gas 
ratio. 

Our CO analyses assumed inner CO fractional abundances 
(with respect to H 2 ) of 1 x 10 “^ for C stars, 6 x 10 "^ for S stars and 
3 X 10 for M stars, in line with canonical results. The sizes of 
the CO envelopes, which are determined by photodissociation, 
were calculated using the results of Mamon et al. (1988 1 . 

In our CO excitation analyses we included radiative tran¬ 
sitions for the hrst 40 rotational energy levels in the ground 


and hrst vibrationally excited states, taken from Chandra et al. 
( |1996| l. The collisional rates, which were only available for the 
ground vibrational state, were taken from Yang et al. ( 2010|l and 


weighted assuming an ortho-/para-H 2 ratio of 3 by Schoier et al. 
(|201 l|l. They cover kinetic temperatures from 2 to 3000 K. 


We calculated the best ht models for each star using a 
statistic, which we dehne as 


2 \ ’ (fmod,/ fobs,/) 

3^ "Zj-5- 


(6) 


where I is the integrated line intensity, cr is the uncertainty in 
the observations (generally assumed to be 20%), and N is the 
number of lines being modelled. We also calculate a reduced 
value such thatx^^d - I~ P) where p is the number of free 
parameters. We t^e the mass-loss rate and the h parameter as 
the two free parameters, i.e. p -2. (Although p is also adjusted, 
it is done so based solely on the line widths and not the complete 
model results.) 

In Eig.[2we show the resultant molecular emission line mod¬ 
els plotted against observational data for two C stars; U Hya, 
which is a low mass-loss rate object, and Al Vol, which is a 
higher mass-loss rate object. Alongside our model results for 
U Hya (shown in blue) we have plotted the results of the same 
model with the only alteration being a constant velocity instead 
of a velocity law following Eq. |^(red dashed line). As can be 
clearly seen, the constant velocity model severely under-predicts 
the higher-7 lines, compared with both the observations and the 
standard model with a velocity prohle. This is a good example of 
both the importance of considering high-7 lines when constrain¬ 
ing a model, and the effect of using a velocity prohle. Eigure]^ 
shows model and observations for Y Lyn, a low mass-loss rate S 
star, and S Cas, a higher mass-loss rate S star. Eig.[^shows model 
and observations for L 2 Pup, a very low mass-loss rate (and low 
expansion velocity) M star, and KU And, a higher mass-loss rate 
M star. 


The remainder of our results are plotted in Eigures A. 1 


A.2 


and IA.3I and the model results are summarised in Table We 
plot new HIEI and IRAM observations where they are available, 
and archival (1 —> 0) and (2^1) lines where new low-7 ob¬ 
servations are missing, to show as broad a range of lines with 
respect to 7 as possible. 


4. Analysis and discussion 

4.1. Modelling target selection 


Not all the stars in our observational sample were modelled in 
this paper. Some stars, such as the semi-regular variables (SRV) 
BK Vir, RX Boo, R Crt, Y UMa, SW Vir, RT Vir, X Her, and 
EP Aqr, will be modelled in a forthcoming paper (Alcolea et al., 
in prep.), while some were excluded because we knew a priori 
that our one-dimensional code, in which spherical symmetry is 
assumed, could not take their idiosyncrasies into account. 

The new observations for unmodelled carbon stars can be 
seen in Pig. B. 1 R Scl was excluded because it is known to have 
a detached shell and a spiral structure in the gas indicative of 
an unseen companion (Maercker et al.|[2012[). See also Schoier 


et al. ( 2005|l for a discussion of difficulties modelling this par¬ 


ticular star. Conversely, U Hya is also known to have a detached 
shell ( Waters et al.|1994] l, but the shell is sufficiently distant that 
the CO in it is most likely photodissociated and hence it has no 
impact on the detected CO lines, enabling us to keep U Hya in 
our modelled sample. IRAS 06192-1-4657 was excluded because 
we had insufficient data for a robust model, particularly since the 
(5 ^ 4) HIEI line was not convincingly detected. TX Psc was 
discounted because it has a two-component molecular wind, and 


is known to have an irregular structure as discussed in Heske 
[etar] ( [T989l l. 
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Table 5. CSE parameters from CO models 


Star 

M 

[Mo yr-‘] 


h 

xii 

N 

C stars 

V701 Cas 

4.5 X 

10 -^ 

0 

5 X 10“ 

5.1 

4 

V384 Per 

2.3 X 

10 -^ 

0 

6 x 10-' 

2.8 

12 

GY Cam 

3.7 X 

10 -^ 

2.0 

9x 10-' 

2.0 

5 

R Lep 

8.7 X 

10 -’ 

2.0 

5 X 10-2 

1.4 

9 

VI259 Ori 

8.8 X 

10 ^^ 

3.0 

5 X 10“ 

1.7 

9 

UU Aur 

1.7 X 

10 -’ 

1.0 

9 X 10-5 

2.3 

8 

V688 Mon 

6.1 X 

10 ^^ 

1.0 

5 X 10“ 

0.58 

5 

Al Vol 

4.9 X 

lO^*^ 

1.0 

9x 10-' 

0.76 

9 

UHya 

8.9 X 

10 -* 

5.0 

8 X 10-2 

2.1 

15 

XTrA 

1.9 X 

10 -’ 

1.0 

5 X 10-' 

1.0 

9 

II Lup 

1.7 X 

10-5 

0 

5 X 10-' 

1.7 

20 

VCrB 

3.3 X 

10 -’ 

1.0 

2 x 10-' 

3.2 

7 

V82I Her 

3.0 X 

10 -^ 

1.0 

7x 10-' 

2.7 

14 

V Aql 

1.3 X 

10 -’ 

0 

5 X 10-' 

2.0 

11 

VI968 Cyg 

7.5 X 

10 -® 

1.0 

8 x 10“ 

2.0 

5 

RV Aqr 

2.3 X 

10-5 

1.0 

5 X 10-' 

3.4 

7 

S stars 

R And 

5.3 X 

10 -’ 

1.5 

6 x 10-' 

1.8 

6 

SCas 

2.8 X 

10-5 

1.0 

5 X 10-' 

2.4 

6 

W And 

2.8 X 

10 -’ 

3.0 

4x 10-' 

0.93 

5 

R Gem 

4.3 X 

10 -’ 

2.0 

6 x 10-' 

2.2 

5 

YLyn 

1.7 X 

10 -’ 

1.5 

5 X 10-2 

0.92 

5 

RCyg 

9.5 X 

10 -’ 

2.0 

1 X 10“ 

1.9 

5 

M stars 

KU And 

9.4 X 

10-5 

2.0 

3 X 10-2 

1.7 

7 

AFGL 292 

2.1 X 

10 -’ 

1.0 

5 X 10-2 

3.5 

4 

RHor 

5.9 X 

10 -’ 

1.0 

4x 10-' 

0.19 

6 

NV Aur 

2.5 X 

10-5 

1.0 

6 x 10-' 

1.8 

11 

BXCam 

4.4 X 

10-5 

1.0 

5 X 10-5 

3.1 

7 

GX Mon 

8.4 X 

10-5 

1.0 

5 X 10-' 

2.6 

14 

Lt Pup 

1.4 X 

10 -® 

0 

2 x 10-' 

1.3 

10 

SCMi 

4.9 X 

10 -* 

0 

6 x 10-' 

1.2 

4 

RLMi 

2.6 X 

10 -’ 

2.0 

2 x 10-' 

2.2 

6 

R Leo 

1.1 X 

10 -’ 

5.0 

5 X 10-2 

5.7 

15 

SCrB 

2.3 X 

10 -’ 

5.0 

2 x 10-' 

3.9 

10 

Vllll Oph 

1.2 X 

10-5 

1.0 

6 x 10-' 

1.9 

10 

RR Aql 

2.4 X 

10-5 

2.0 

5 X 10-2 

2.7 

7 

VI943 Sgr 

9.9 X 

10 -® 

1.0 

5 X 10-5 

3.1 

4 

VI300 Aql 

l.Ox 

10-5 

1.0 

6 x 10-2 

4.6 

16 

TCep 

9.1 X 

10 -® 

1.0 

3 X 10-' 

0.23 

3 


Two S stars were excluded from modelling, RS Cnc and 
Gru, both of which have two velocity components in their line 
profiles, as can be seen in Fig. B.2 RS Cnc is known to have a 
bipola r outflow and a disc structure around the star ( [Libert et al. 


2010 1 , and Gru is known to have a bipolar outflow, a GOV 


binary companion, and evidence of a second hidden companion, 
as described by Mayer et al. (2014|l. 


The unmodelled M-type star sample contains the SRVs listed 
above and R Hya and V370 And, which were excluded for hav¬ 
ing double component winds. 

These exclusions have left us with a sample for which we 
have plotted histograms of mass-loss rates and expansion ve¬ 
locities by chemical type in Fig. Note in particular the gap 
in expansion velocity for M stars with no modelled stars in the 
range 10 kms^^ < i/oo < 15 kms *. There is a similar gap in the 
S star distribution, however this is more a result of having only 
one S star of high mass-loss rate and high velocity. The C stars 
are fairly evenly distributed in velocity. 







Velocity [km/s] Velocity [km/s] 



Velocity [km/sl 


Fig. 1. Models (blue) and observed data (black) of two exam¬ 
ple C stars, plotted with respect to LSR velocity. An * next to 
the telescope name indicates that archival data is plotted. The 
dashed red lines on U Hya show the predicted model if a con¬ 
stant expansion velocity is used, keeping all parameters the same 
as in the blue model. 


4.2. Goodness of fit 


The goodness of fit of our models has been primarily estimated 
using the;^^^ method discussed in Sect. 3.2 To further visualise 
the goodness of fit we have plotted the values of modelled line 
intensity divided by the observed line intensity (/mod/fobs) for 
each observed line in each star. The resulting plots can be seen 
in Fig. 


A.4 


for C stars. Fig. A.5 for S stars and Fig. A.6 


forM 

stars. For the large majority of sources, we can model the CO 
lines with upper energy levels below 250 K (and below 580 K 
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Fig. 2. Models (blue) and observed data (black) of two example 
S stars, plotted with respect to LSR velocity. 


for sources with observed J — 14 —» 13 lines) with reasonably 
accuracy using a constant mass-loss rate. In particular, there is 
no trend with 7-number. 

We have also combined the same quantity for all stars 
grouped by line, for the (1 ^ 0), (2 ^ 1), (5 ^ 4), and (9 —> 8) 
lines. The resultant plots can be seen in Fig.[^ The distributions 
are reasonably symmetric, except for the (1 —> 0) that seems to 
be under-predicted in the models for all three chemical types. 
The (1 —» 0) line intensity is particularly sensitive to the size 
of the CO envelope and the former discrepancy can be (at least 
partly) remedied if a larger CO envelope is used than the size 
predicted by the Mamon et al. ( |1988| l model. 

However, this line, and to some extent the (2 ^ 1) line 
also have a tendency to towards double-peaked prohles rather 
than the flat or slightly rounded observed line prohles. This hap¬ 
pens in about 5 out of 38 cases and is a well-known fact in CO 
line modelling, for example see |01ofsson et al. ( |2002| l. An in¬ 
crease in the size of the CO envelope would further enhance the 
double-peaked nature of the line prohles. In addition to the CO 
envelope size, there are several possible reasons for the double- 
peaked (1^0) line prohles: a too-warm outer CSE, a too- 
distant source, or maser action in the inner CSE (which was in 
evidence in some of our models but is known to be produced in 




Velocity [km/s] Velocity [km/s] 


Fig. 3. Models (blue) and observed data (black) of two example 
M stars, plotted with respect to LSR velocity. An * next to the 
telescope name indicates that archival data is plotted. 


nature). Since we have studied a large sample of sources in this 
paper, we have preferred to treat all sources the same and, there¬ 
fore, make no adjustments for individual sources, consequently, 
we conclude that we produce overall good hts to the CO line 
intensities for the majority of our sources, while there remain 
some discrepancies (some known from previous studies, such as 
the tendency for double-peaked model line prohles), we believe 
they have no major impact on the conclusions in this paper. 

4.3. Trends with mass-loss rate and chemical type 

To investigate any trends that appear in our results, we plot¬ 
ted mass-loss rate against luminosity, expansion velocity, opti- 
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10 
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Fig. 4. The numbers of stars with the indicated expansion veloc¬ 
ities and mass-loss rates. C stars are grey, S stars are white and 
M stars are black. 



differences in wind-driving efficiency between the two chemical 
types. 

The plot of mass-loss rate against h parameter does not show 
any obvious trends, indicating that the dust properties embedded 
in the h parameter do not directly depend on the mass-loss rate. 

We have modelled CO lines covering a relatively broad range 
in energies, which means that the kinetic temperature distribu¬ 
tion is well constrained, except perhaps for the very inner part. 
In Fig. 1^ (bottom right) we plot the kinetic temperature in the 
CSE at a radius of 100 x against mass-loss rate. There is no 
obvious correlation, but the trend is that the CSEs become, on 
average, cooler the higher the mass-loss rate. This is certainly 
expected, since, with all else being equal, a high mass-loss rate 
leads to more efficient CO line cooling and a lower drift velocity 
contributing to less efficient heating of the gas. 


4.4. Comments on individual stars 


Eor the discussion of individual stars below, note that the mod¬ 
elled and observed lines plotted in Eigures[T and|^ and in the 
appendices in Eigures |A.l A.2[ and |A.3[ show a maximum of 
five emission lines. In most cases there are several archival ob¬ 
servations which were included in the modelling, and are listed 
in Table C. 1 but are not plotted with the new observations. These 
archival observations can have a significant impact on the choice 
of best fit model. The full collection of line observations for each 


star is included in the goodness of fit plots in Eigures A.4 A.5 


and IA.6I The number of lines included in the radiative transfer 
modelling for each star is noted in Table 


4.4.1. C stars 


Fig. 5. Goodness of fit as defined by model/observed integrated 
line intensity of all stars per transition line for C (grey), S (white) 
and M (black) stars. N is the number of lines of each chemical 
type in each bin. 


cal depth at 10 fim, and the h parameter in Fig. There is a 
clear trend between mass-loss rate and luminosity, which does 
not seem to depend on chemical type. This is in line with the 
expectation that the higher luminosity stars should have higher 
mass-loss rates, either due to having a higher mass, or due to 
a more advanced age on the AGB. There is also a correlation 
between expansion velocity and stellar luminosity, although it is 
less tight than the mass-loss rate and luminosity correlation. This 
suggests that the gas expansion velocity is tied to the stellar lumi¬ 
nosity. It is not surprising, then to find a correlation between the 
mass-loss rate and gas expansion velocity. The relation between 
these two mass-loss characteristics puts constraints on any vi¬ 
able mass-loss mechanism. Although there is some segregation 
in both of these plots, this is due to gaps in our sample rather 
than intrinsic trends. 

There is also a correlation between mass-loss rate and opti¬ 
cal depth at 10 fim, as is expected if the gas and dust mass-loss 
rates are correlated. The correlation has slightly different slopes 
for the C and M stars. For C stars we find the best fit slope is 
1.0 + 0.3 and for M stars 1.3 + 0.6. There are not enough S stars 
in our sample for a meaningful determination. It should also 
be noted that the optical depth is taken at 10 /rm, near a strong 
silicate feature in the M stars and an SiC feature in the C stars. 
This probably contributes to the difference in slope, along with 


VI259 Ori, V688 Mon, AI Vol, V821 Her, VI968 Cyg, and RV 
Aqr all have higher mass-loss rates (above 10“® Mq yr“') and 
expansion velocities in excess of 10 kms“*. Their CO line pro¬ 
files can all be described by models that include a velocity pro¬ 
file as described in Eq. with jS = 1, except for in the case of 
V1259 Ori for which /3 = 3. V688 Mon has ISM contamination 
in the low-7 data for which we compensated when calculating 
the integrated line intensities. 

UU Aur, U Hya, X TrA, and V CrB are the low mass-loss rate 
carbon stars, all having mass-loss rates lower than 10 “® Mq yr“' 
and expansion velocities <10 km s *. Their CO line profiles are 
all reasonably well-described by our models with /I = 1, except 
for in the case of U Hya for which /3 - 5. C Hya also has the 
lowest mass-loss rate (8.9 x 10“^ Mq yr“^) of the carbon stars in 
this sample. 

There were a handful of stars for which a constant velocity 
model (J3 - 0) gave a markedly better fit than the standard ve¬ 
locity profile model. The stars which fell into this group were 
V701 Cas, V384 Per, V Aql, and II Lup. The first two objects 
have moderate mass-loss rates and V Aql has a relatively low 
mass-loss rate for a C star. The models of these three objects 
all suffer from maser emission in the low-7 lines if a velocity 
profile is used (which does not entirely disappear in V Aql even 
with a constant velocity as can be seen in Fig. |A.l| i. The model 
maser emission arises in the inner regions of the CSE. It may 
be that, for various reasons, such maser action is not produced 
in nature. V701 Cas also suffers from ISM contamination in the 
low-7 lines which was compensated for. II Lup is both a high 
mass-loss rate and high expansion velocity object. The line in¬ 
tensities of a jS > 0 model do not match the observed data as well 
as a constant velocity model by a significant amount: the best fit 
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Fig. 6. Mass-loss rate and expansion velocity plotted against luminosity; mass-loss rate plotted against expansion velocity, optical 
depth and the /i-parameter; temperature in the CSE at a radius of 100 x plotted against mass-loss rate. Carbon stars are red 
circles, S stars are green squares, M stars are blue triangles. 


13-1 model has = 3.3 compared with x^^^ = 1.7 for the 
constant velocity model. 

GY Cam has a narrow component (~ 8 kms“^ wide) lying 
approximately at the systemic velocity of i/isr = -49 kms“* that 
is clearly visible in the (1 —» 0), (2 —» 1) and (5 —» 4) lines. Note 


that it is absent in the (9 —» 8) line, which suggests that it may 
be cooler in temperature than the gas contributing to the broader 
component of the line profile. Our CO model was not able to 
take the narrow component in these lines into account and we 
did not compensate for it in any way. Since our best-fit models 
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are determined from total line intensity, this had a slight effect 
on our result. However, it does not seem to have been as large an 
issue as in some stars such as AFGL 292 (see Sect.|4.4.3lbelow). 


4.4.2. S stars 


The majority of the S stars were straightforward to model. 
Although R And, R Gem, and Y Lyn are known to have binary 
companions ( Proust et al.| 19811 |Pourbaix et al.||2003| l, this did 
not seem to have any significant impact on the observed line pro¬ 
files when compared to the line models. The most problematic 
of these five stars was R Gem, which has the noisiest data. 


The only S star which we found problematic to model was 
S Cas. It has the largest terminal expansion velocity (more than a 
factor of two higher than any of the other S stars) and we found it 
to also have the highest mass-loss rate. It was also the only S star 
for which we detected any CN “bonus” lines (see Sect. |B.2| . In 
these ways it bears many similarities to some of the carbon stars. 


4.4.3. M stars 


KU And, NV Aur, BX Cam, GX Mon, Vllll Oph, and 
V1300 Aql are the oxygen-rich stars in our sample with the high¬ 
est terminal expansion velocities (Voc > 15 kms^'). They can all 
be described well by models which include a velocity profile 
with jS = 1 except for NV Aur, for which /3 - 2. It should also 
be noted that the NV Aur lines exhibit some ISM contamination, 
especially in the low-7 lines. 

R Hor, R LMi, S CrB, RR Aql, V1943 Sgr, and T Cep 
have relatively low terminal expansion velocities, with Voc < 
10 kms *. Their CO profiles are described well by our models 
and have velocity profiles with p in the range 1.0</3<2.0. 

L 2 Pup and S CMi both appear to have very low terminal ex¬ 
pansion velocities, with Voc - 2 kms *, and we have therefore 
modelled the CO line emission in both cases assuming a con¬ 
stant velocity. The fits to the observed lines are in general very 
good, although it is unclear whether the S CMi (1^0) line 
is really as strongly double-peaked as it appears, in which case 
this behaviour is not reproduced in our model. However, there 
are indications that L 2 Pup has a central toroidal or disc struc¬ 
ture (Lykou et al. 2015[ [Kervella et'aL]|2014|l, with a possible 
close companion and bipolar outflow ( [Kervella et al.||2015] l. In 
this case, our spherically symmetric model most likely does not 
represent a realistic view of L 2 Pups circumstellar environment, 
and the narrow line is not produced in a slowly, spherically ex¬ 
panding envelope. Similarly, Ragland et al. ( |2006| l have found 
evidence of asymmetry in S CMi, suggesting non-spherical ge¬ 
ometry. 

AFGL 292 appears to have a two-component line profile. 
The new IRAM observations in particular suggest an additional 
narrow component centred at the stellar systemic velocity. This 
component is not visible in the high-7 lines observed with HIFI, 
which could be an effect of cooler gas or simply because the 
signal-to-noise ratio of the HIFI observations is insufficient. As 
a result, the model we present in Fig. A. 3 is the best fit based 
on the integrated intensities, but does not match the line profiles, 
particularly not for the low-7 lines. 

R Leo shows asymmetric line profiles, with the HIFI lines 
being less bright on the red-shifted side, as can be seen in Fig. 
|A.3| Conversely, some of the low-7 lines, most notably the (2 —> 
1) IRAM line and the (3 ^ 2) CSO line, have significant peaks 
on the red side of the line profiles. Also, the (1 —> 0) IRAM 
line seems to have a central peak as well as two outer peaks. 



M comparison: Other studies / This study 


Fig. 7. Comparisons with mass-loss rates from past studies: 
■^previous/-^new. grouped by chemical type. C stars are grey, M 
stars are black and S stars are white. 


This could be indicative of localised clumps of gas at different 
temperatures. 


4.5. Comparison with other studies 


We compared our modelling results with those of other studies 
of stars in our sample. Correcting the mass-loss rates for both 
different distances to the object and different CO abundance as¬ 
sumptions, we compared results using the metric Mprevious/47new, 
where “previous” indicates the mass-loss rate found in an earlier 
study and “new” refers to our mass-loss rate results. The previ¬ 
ous studies with the largest overlap with our sample were those 


ot'ISchoier & Olot'sson (2001|l, Knapp et al.|(1998|l, De Beck 

et al.|(2012||, and Ramstedt & Olofsson|(20141l. 

Where possible 
ir, coming from 

we also included an extra data point for each sta 

other overlapping studies such as Olofsson et al. (1993 

,| Woods 

et al.|(20031l, Neri et al. ( 1998|l,|Ramstedt et al.|( 

2008b, 

Teyssier 


et al.| ( |20()^ , |Young| ( |1995| l, md |01ofsson et al.| ( |2002| l. Grouping 


all the comparison ratios into chemical type, we find the mean 
ratio and standard deviation for C stars is 1.4 + 0.9, for M stars 
is 1.4 ± 0.8 and for S stars is 1.9 ± 0.8. The collected ratios are 
also plotted in histograms in Fig.|^ 


In general, our models tend to give mass-loss rates on aver¬ 
age 40% lower than previous studies, even after correcting for 
distances and and CO abundance assumptions. Apart from this 
offset, the spread in the estimates is about +50%. We present 
the first study where high-7 lines up to (9^8) — or in some 
cases (14 —» 13) — are used in the analysis of a large sample 
of stars. It could hence be inferred that taking high-7 lines prop¬ 
erly into account when modelling mass-loss rates results in lower 
predicted mass-loss rates. Another possible factor is in our use 


and shown in Fig.[]J the inclusion of an accelerating wind model 
compared with a constant expansion velocity wind can have a 
significant effect on the resultant line intensities (and intensity 
ratios between different lines). 


of the velocity profile described in Eq.H As discussed in 3.2.2 
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5. Conclusions 

We have presented new HerschellHWl and IRAM CO line data 
for a sample of 53 AGB stars, which cover a wide range of tran¬ 
sitions from J - 1 —>0 up to 7 = 9 —>8 (and J - 14 —> 13 in 
a few cases). For 38 of the stars, we used radiative transfer mod¬ 
els to determine their mass-loss rates. Our results cover a wide 
range of mass-loss rates, from ~ - 2 x 10 “^ Mq yr“', and 

include all three chemical types of AGB stars (M, S, and C). In 
general, we find it possible to get a very good fit to the observed 
CO lines which cover up to 250 K in energy (in a few cases 
up to 600 K) without the need to invoke mass-loss rate changes 
with time. Our model results are also in reasonable agreement 
with past studies (within a factor two for 63% of the stars), al¬ 
though they generally predict slightly lower mass-loss rates than 
past studies based primarily on lower-7 lines. We found that our 
models under-predict the CO (1 —> 0) line intensity more often 
than other lines, although the precise reason for this is not clear. 

In analysing our results we found correlations between mass- 
loss rate and stellar luminosity, gas terminal expansion veloc¬ 
ity, and dust optical depth. The latter correlation seems to show 
slightly different slopes for the M and C stars, with the S stars 
falling in between the two. We also found that there is a tendency 
for the CSEs to have a lower kinetic temperature the higher the 
mass-loss rate, an effect expected due to less efficient dust-gas 
collision heating and stronger CO line cooling in the denser ob¬ 
jects. 

We also found that the inclusion of an expansion velocity 
profile in our models can have a significant effect on the results. 
In particular, the use of a constant velocity instead of an acceler¬ 
ating wind can severely under-predict the high-7 line intensities, 
such as those observed by HerschelfilWl. The exact shape of the 
velocity profile is obtained by fitting the shapes of, in particular, 
the higher-7 lines. 

Our observations also included a series of “bonus” lines, 
which fell within the observing ranges of our CO observations. 
'^CO was the most commonly detected “bonus” line in all three 
chemical types of AGB stars. It was detected in eight C stars, 
five S stars and eight M stars. CN and HC 3 N were detected in 
C stars, with one S star (S Cas) also being detected in some of 
the observed CN lines. SiO was frequently detected in M stars 
(17 detections), as opposed to the C stars where only two were 
detected. SiS was detected in the higher mass-loss rate C and 
M stars, although not in all of the highest mass-loss rate stars. 
There were seven detections in C stars and five in M stars. 
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Table 2. HIFI and IRAM CO line observations 


Star 

Voo 

[km s“^] 

IRAM 

CO (1 ^ 0) CO (2 ^ 1) 

[Kkms“‘] [Kkms“‘] 

CO (5 ^ 4) 

[K km s“'] 

HIFI 

CO (9 ^ 8 ) 

[K km s“*] 

CO (14^ 13) 
[K km s“'] 

C stars 

RScl 

16.5 

81.9 (1.0) 

103.3 (0.9) 

5.30 (0.08) 

1.57 

(0.07) 

1.33 (0.38) 

V701 Cas 

11.5 

*36.4 (0.3) 

* 44.4 (0.4) 

1.67 (0.22) 

1.93 

(0.07) 


V384 Per 

15 

63.6 (0.2) 

90.1 (0.3) 

5.52 (0.18) 

6.25 

(0.09) 


GY Cam 

20 

43.8 (0.2) 

70.1 (0.3) 

3.14 (0.20) 

2.90 

(0.08) 


R Lep 

18 

32.5 (0.4) 

100.3 (0.6) 

4.23 (0.13) 

4.57 

(0.08) 


VI259 Ori 

16 

63.8 (0.2) 

85.7 (0.3) 

3.93 (0.21) 

3.23 

(0.07) 


IRAS 06192+4657 

6 



< 0.22 (0.16) 

0.290 

(0.171) 


UU Aur 

12 



2.47 (0.15) 

2.31 

(0.08) 


V 688 Mon 

13.5 

*30.8 (0.3) 

*67.4 (0.2) 

1.96 (0.19) 

2.31 

(0.09) 


AI Vol 

12 



6.87 (0.13) 

8.12 

(0.08) 

7.62 (0.23) 

UHya 

6.5 

8.87 (0.37) 

53.7 (1.0) 

2.98 (0.16) 

2.07 

( 0 . 11 ) 


XTrA 

6.5 



1.36 (0.18) 

1.76 

(0.08) 


II Lup 

21.5 



16.8 ( 0 . 1 ) 

17.6 

(0.07) 


VCrB 

7.5 



0.541 (0.159) 

1.07 

(0.08) 


V821 Her 

13.5 

82.8 (0.5) 

136.5 (0.7) 

6.32 (0.15) 

7.46 

(0.08) 


V Aql 

11 

12.6 ( 0 . 2 ) 

23.1 (0.3) 

1.36 (0.13) 

1.30 

(0.07) 


VI968 Cyg 

20 

37.2 (0.3) 

1 

5.58 (0.18) 

5.79 

(0.07) 


RV Aqr 

15 

46.0 (0.2) 

73.9 (0.4) 

3.64 (0.13) 

4.16 

(0.07) 


TXPsc 

4 



0.887 (0.141) 

1.24 

(0.08) 


S stars 

R And 

8 

27.0 (0.4) 

56.9 (0.4) 

3.72 (0.16) 

5.58 

(0.08) 


S Cas 

19 

37.6 (0.3) 

90.9 (0.5) 

4.27 (0.18) 

6.08 

(0.08) 


W And 

6 

8.14 (0.1) 

24.5 (0.4) 

1.07 (0.13) 

1.34 

(0.09) 


R Gem 

5 

7.01 (0.3) 

10.9 (0.6) 

0.565 (0.230) 

0.471 

(0.156) 


YLyn 

8 

11.1 (0.4) 

29.6 (0.6) 

1.52 (0.19) 

1.28 

(0.09) 


RS Cnc 

2.5 

22.9 (0.2) 

63.4 (0.5) 

3.60 (0.19) 

3.17 

(0.08) 


RCyg 

9 

15.6 (0.3) 

43.2 (0.5) 

2.21 (0.16) 

1.42 

(0.08) 


7t Gru 

10 



13.8 (0.15) 

13.9 

(0.04) 

10.0 ( 0 . 1 ) 

M stars 

KU And 

20 

48.9 (0.3) 

83.5 (0.4) 

3.15 (0.16) 

4.41 

(0.07) 


V370 And 

9 



2.96 (0.22) 

2.35 

(0.09) 


AFGL 292 

8.5 

4.56 (0.23) 

10.1 ( 0 . 2 ) 

0.492 (0.229) 

0.481 

(0.127) 


RHor 

4 



4.23 (0.14) 

3.31 

(0.08) 


NV Aur 

18 

43.3 (0.2) 

59.1 (0.2) 

2.54 (0.13) 

1.84 

(0.07) 


BX Cam 

19 

33.9 (0.2) 

50.8 (0.2) 

2.96 (0.19) 

2.58 

(0.08) 


GXMon 

19 

64.2 (0.1) 

128.1 ( 0 . 2 ) 

5.83 (0.18) 

4.07 

(0.07) 


L 2 Pup 

2 



1.82 (0.14) 

2.56 

(0.09) 

3.41 (0.24) 

S CMi 

2 

0.608 (0.358) 

3.84 (0.60) 

0.271 (0.108) 

0.436 

(0.084) 


RLMi 

7.5 



1.51 (0.13) 

1.27 

(0.08) 


R Leo 

8.5 



5.86 (0.15) 

8.35 

(0.08) 


RCrt 

12 



3.57 (0.18) 

3.26 

(0.08) 


BK Vir 

6 

3.73 (0.16) 

13.1 (0.2) 

0.474 (0.101) 

1.03 

(0.08) 


YUMa 

7.5 



0.735 (0.088) 

0.807 

(0.090) 


RT Vir 

8.5 



1.90 (0.14) 

1.53 

( 0 . 10 ) 


SW Vir 

8.5 



4.34 (0.23) 

4.14 

(0.09) 


RHya 

10 



6.94 (0.22) 

7.85 

(0.09) 

3.98 (0.17) 

RXBoo 

10 



4.52 (0.23) 

3.79 

(0.09) 

2.92 (0.24) 

SCrB 

7 

3.19 (0.33) 

! 

1.19 (0.11) 

0.910 

(0.092) 


XHer 

8.5 



2.37 (0.18) 

1.98 

(0.08) 


Vllll Oph 

17 

47.2 (0.3) 

82.1 (0.4) 

3.13 (0.19) 

2.50 

(0.07) 


RR Aql 

9 



1.50 (0.11) 

1.14 

(0.08) 


VI943 Sgr 

6.5 



0.862 (0.168) 

0.754 

(0.082) 


VI300 Aql 

18 

68.9 (0.4) 

102.7 (0.7) 

3.66 (0.15) 

3.12 

(0.06) 


TCep 

5.5 



1.44 (0.10) 

1.90 

(0.09) 


EP Aqr 

12 



3.73 (0.16) 

3.35 

(0.08) 



Notes. The value in brackets after the flux gives the integrated noise RMS. An ellipsis (...) indicates that the line was not observed for the indicated 
star; * indicates that flux has been corrected for ISM emission; ! indicates the line was observed but with unreliable flux calibration. 
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Table 4. SED parameters 


Star 

D 

L. 

T. 

fio 

Reference 


[pc] 

[Eol 

[K] 


[xlO*^ cm] 


C stars 


V701 Cas 

1720 

7800 

2800 

0.85 

2.1 

Menzies et al. 

20061 


V384 Per 

560 

8100 

2000 

0.25 

2.1 

t 

Schoier et al. 

i2U13 

i 

GY Cam 

1030 

7800 

2000 

0.3 

2.0 

G 

roenewegen et al. 

(2 

002b 

R Lep 

432 

5500 

2200 

0.06 

1.7 

T 

Ramstedt & Olotsson 

^014 

V1259 0ri 

1600 

9300 

2200 

1.7 

2.2 

IV 

lenzies et al. 

TDDSr 


UU Aur 

260 

6900 

2800 

0.017 

1.9 

T 

Schoier et al. 

iTOTBi 

V688 Mon 

1770 

8800 

2800 

0.6 

2.2 

Menzies et al. i 

2DD6r 


AI Vol 

710 

9000 

2100 

0.45 

2.2 

TISchoier et al. 

1 2007 

( 

UHya 

208 

*4000 

2400 

0.012 

1.5 

tiRamstedt & Olotsson 

(2014 

XTrA 

360 

5400 

2200 

0.024 

1.7 

Cox et al. 


Zi 



II Lup 

500 

8800 

2400 

0.55 

2.2 

t 

Schoier et al. 

(20131 

VCrB 

630 

5300 

1800 

0.035 

1.7 

Cox et al. 

liUL 

U 



V82I Her 

600 

7900 

2200 

0.45 

2.1 

T 

Schoier et al. 

i2013 

1 

V Aql 

330 

6500 

2800 

0.02 

1.9 

t 

Schoier et al. 

(2013 

i 

VI968 Cyg 

1480 

10200 

2400 

0.85 

2.3 

IV 

lenzies et al. i 

2DD^ 


RV Aqr 

670 

6800 

2200 

0.27 

1.9 

t 

Schoier et al. 

(2013 

( 


S stars 


R And 

350 

6300 

1900 

0.05 

1.8 

t 

Ramstedt & Olofsson 

2014 

S Cas 

570 

8000 

1800 

0.5 

2.1 

t 

Ramstedt & Olotsson 

2014 

W And 

450 

5800 

2400 

0.1 

1.8 

van Ceeuwen 

12DD7 

( 



R Gem 

820 

5500 

2400 

0.035 

1.7 

Whitelock et 

al. 

m 

08 


YLyn 

253 

*4000 

2400 

0.02 

1.5 

Ramstedt & Ulo 

tsson 

2014 

RCyg 

690 

6200 

2600 

0.14 

1.8 

Whitelock et al. 

^2008 

1 


M stars 


KU And 

680 

11800 

2000 

0.90 

2.5 

t 

Schoier et al 

(2013( 

AFGL 292 

319 

*6000 

2200 

0.02 

1.8 

Winters et al. ( 

2003 

RHor 

310 

8500 

2200 

0.30 

2.1 

T 

Schoier et al 

ITOTSi 

NV Aur 

1200 

9800 

2000 

3.50 

2.3 

t 

Schoier et al 

(2013 I 

BXCam 

500 

7500 

2800 

1.30 

2.0 

t 

Schoier et al 

(2013 I 

GXMon 

550 

8200 

2600 

2.00 

2.1 

t 

Schoier et al 

(2013 I 

Lt Pup 

86 

*4000 

2800 

0.07 

1.5 

Schoier et at. i 

ion I 

S CMi 

470 

5000 

2800 

0.07 

1.6 

Knapp et at. 1 1 


' 

RLMi 

330 

5500 

2400 

0.2 

1.7 

Whitetock et a 

1 . 1 

D08 

R Leo 

130 

4600 

1800 

0.10 

1.6 

I 

Schoier et al 


SCrB 

400 

5400 

2400 

0.2 

1.7 

Knapp et al. (|l 

^1^8 

Vllll Oph 

750 

7500 

2000 

0.75 

2.0 

T 

Schoier et al 

fimsi 

RR Aql 

530 

7900 

2000 

0.70 

2.1 

t 

Schoier et al 

(pUTT( 

VI943 Sgr 

200 

5000 

2200 

0.05 

1.6 

van Leeuwen ( 

2007: 

VI300 Aql 

620 

10600 

2000 

3.50 

2.4 

T 

Schoier et al 

(12001 

TCep 

190 

5700 

2400 

0.10 

1.8 

van Leeuwen ( 

20071 


Notes. An * indicates that the luminosity is assumed (for semi-regular variables) rather than calculated. The references listed are for distances, 
D, and in some cases the luminosity, L, where no period is known. A t in the references indicates that all parameters including the effective 
temperature, T„ and optical depth, Tio were also taken from the referenced DUSTY modelling. The absence of a t indicates that temperature and 
optical depth were calculated as part of this work.indicates that the distance was calculated from the period-magnitude relation in the cited work. 
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Appendix A: Modelled stars 

Here we present the observations and models for the stars not 
included in the body of the paper. In each instance, we present 
all new data from IRAM and HIFI — black histograms — over¬ 
plotted with model results for the parameters given in Table 
For those stars that were not observed with IRAM, we include 
archival CO (1 —» 0) and (2 —» 1) lines from various telescopes 
as available. These archival lines are indicated by an * next to the 
telescope name in the plot and allow us to present an overview of 
our models from low- to high-/. R Hor is the only star for which 
these low-T lines were not available. Our model for R Hor still 
incorporates some low-/ lines as noted in Table [CT] 

The plots for C stars are shown in Fig. A. 1 for S stars in Fig. 


A.2 |and for M stars in Fig. |A.3| 

As well as plotting the model and observed CO lines, we 
have calculated “goodness of ht” per line in each star. This gives 
us an indication of which lines may be outliers or whether there 
are any trends across lines. 

Goodness of fit for C stars is shown in Fig. A.4 for S stars in 
Fig. A.5 and for M stars in Fig. A.6 See Fig.and Sect. 4.2 for 


a discussion of goodness of ht across the entire sample. A list of 
archival lines included in our models is given in Table |C.l1 
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Fig.A.l. Models (blue) and observed data (black) of C stars, plotted with respect to LSR velocity. An * next to the telescope name 
indicates that archival data is plotted. 
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Fig. A.l. cont. Models (blue) and observed data (black) of C stars, plotted with respect to LSR velocity. An * next to the telescope 
name indicates that archival data is plotted. 
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Fig. A.l. cont. Models (blue) and observed data (black) of C stars, plotted with respect to LSR velocity. An * next to the telescope 
name indicates that archival data is plotted. 
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Fig. A.2. Models (blue) and observed data (black) of S stars, plotted with respect to LSR velocity. 
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Fig. A.3. Models (blue) and observed data (black) of M stars, plotted with respect to LSR velocity. An * next to the telescope name 
indicates archival data is plotted. 
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Fig. A.3. cont. Models (blue) and observed data (black) of M stars, plotted with respect to LSR velocity. An * next to the telescope 
name indicates archival data is plotted. 
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Fig. A.3. cont. Models (blue) and observed data (black) of M stars, plotted with respect to LSR velocity. An * next to the telescope 
name indicates archival data is plotted. 
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Fig.A.4. Goodness of fit as defined by model/observed intensity for C stars. 
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Fig. A.5. Goodness of fit as defined by model/observed intensity for S stars. 
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Fig. A.6. Goodness of fit as defined by model/observed intensity for M stars. 
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Appendix B: New observations 

B.1. CO lines 

Of the data described in Table those stars for which we ran 
radiative transfer models were plotted in the body of the paper 
and in Appendix The remaining lines, which were excluded 
from modelling for various reasons (see discussion in Sect. 4.1| l 
are now presented here. The C stars are plotted in Fig. B.l the 
S stars are plotted in Fig. B.2 and the M stars are plotted in Fig. 


In particular, the unusual line profile due to the presence of 
a detached shell can be seen in the C star R Scl and double¬ 
component winds are clearly evident in C star TX Psc and S 
stars RS Cnc and Gru. 

The observation identifiers (ObsIDs) for our Herschel obser¬ 
vations are listed in Table lB.il 


B.2. Bonus lines 


As mentioned in Sect. 2.4 we acquired some “bonus” line spec¬ 


tra for molecules that were observable within our target fre¬ 
quency ranges. In HIFI, in the same range as the CO (5 ^ 4) 
line, we detected SiO (13 —> 12) at 564.249 GHz. As can be 
seen in Table |B.2| it was mostly detected in M stars, especially 
those of lower mass-loss rates, which is in agreement with the 
trend found by [Gonzalez Delgado et al. (2003 1 and the calcula¬ 
tions performed by Schoier et al. 1 2004|l. 


Our detections are plotted in Fig. B.4 There was one detec¬ 
tion in an S star, RS Cnc, which is the most “M-like” S star in 
our sample, based on optical classifications. There were also two 
detections in C stars: V384 Per and V821 Her. They both have 
mass-loss rates in the range ~ 2-3 xl0“® Mq yr“*, putting them 
in the mid-to-high mass-loss rate range. They are located at 560 
and 600 pc respectively, making them two of the nearest C stars 
in the higher mass-loss rate range (> 10“® Mq yr“^). This could 
be why they had (weak) detections, while there were no detec¬ 
tions in other C stars. The two C stars are among the sample 
Schoier et al.j ( (2006] l. These authors used 
ines from 7 = 8 —> 7 down to J - 2 ^ 1 


modelled in SiO by 
observations of SiO 


and there are six overlapping stars between their sample and the 
one in this paper, leaving four stars (Al Vol, II Lup, RV Aqr, and 
R Lep) detected in the lower-7 SiO lines but not in the higher-7 
HIFI line. 

Covered by our IRAM observations was the '^CO (1 —> 0) 
line at 110.201 GHz. The integrated intensities for these detec¬ 
tions are given in Table B.2 and the observations are plotted in 
Fig. 


B.5 


The ‘^CO (1^0) line seems to have been most re¬ 
liably detected in higher mass-loss rate sources across the three 
chemical types. It was not detected at all in stars with mass-loss 
rates ~ 10“®-10“^ Mq yr"^ (note, however, that we are only 
dealing with two stars in this range) and was detected increas¬ 
ingly often for increasing mass-loss rates across chemical types. 

The SiS (6 5) line at 108.924 GHz was also detected in 

seven C stars, five M stars and no S stars. The integrated inten¬ 


sities are listed in Table B.2 and the spectra are plotted in B.6 


The C stars with detections were all in the mass-loss rate range 
~ 10“® to 10“^ Mq yr“', with no detections for lower mass- 
loss rate objects and only one detection out of the two highest 
mass-loss rate C stars observed. Of the M stars, SiS was also de¬ 
tected in the higher mass-loss rate objects, but not in the highest 
mass-loss rate star, V1111 Oph. This trend suggests that SiS is 
more readily formed — or at least more readily detectable — in 


sources of intermediate mass-loss rate, around the range ~ 10 ° 
to 10“^ Moyr^'. 

CN lines were covered by both HIFI and IRAM observa¬ 
tions. The CN (59/2 ^ 47 / 2 ) and (5ii/2 ^ 49 / 2 ) line groups with 
rest frequencies taken as 566.693 GHz and 566.947 GHz were 
covered in the observing range for CO (5 —> 4) and were de¬ 
tected in a handful of C stars. Our IRAM observations covered 
the CN A = 1 ^ 0 lines at 113.123 GHz and 113.488 GHz for 
the (I 1 / 2 -O 1 / 2 ) and (I 3 / 2 -O 1 / 2 ) line groups, respectively, and the 
CN A = 2 ^ 1 lines at 226.617 GHz, 226.874 GHz and 226.360 
GHz for the ( 23/2 - I 1 / 2 ), ( 23/2 - I 3 / 2 ), and ( 25/2 - I 3 / 2 ) line 
groups, respectively. The hyperfine structure of the low-A CN 
lines can be seen particularly clearly. The integrated intensities 
of each line group are given in Table B.3 and the observations 


themselves are plotted in Fig. B.7 


Low-A CN lines were detected in all of the observed C stars. 
Not all lines were detected in all stars, however. The lowest 
mass-loss rate star, U Hya, did not yield a clear detection of 
the (I 1/2 - O 1 / 2 ) or ( 23/2 - I 3 / 2 ) groups, although the remain¬ 
ing lines, including the ( 69/2 ^ 47 / 2 ) and (5 11/2 —> 49 / 2 ) groups 
were clearly seen. The ( 2 ^/ 2 — 13 / 2 ) was also not detected in V701 
Cas, V1259 Ori, V688 Mon, or V821 Her, all of which are rela¬ 
tively high mass-loss rate objects with M ~ 10“® -10“^ Mq yr"^. 
One S star, S Cas, was also detected in CN, in the (I 3/2 - O 1 / 2 ), 
( 23/2 - I 1 / 2 ), ( 23/2 - I 3 / 2 ), and ( 25/2 - I 3 / 2 ) line groups. S Cas 
is the highest mass-loss rate and expansion velocity S star and, 
from its optical classification of S4/6e, is on the higher C/O end 
of the S star scale. 

The last bonus line we detected was HC 3 N (12 —> 11) at 
109.174 GHz. The integrated intensities for the detections are 
listed in Table B.3 and the spectra are plotted in Fig. B.6 HC 3 N 
was only detected in C stars and not in the three lowest mass- 
loss rate objects with mass-loss rates below 10“® Mq yr“^ This 
is probably due to a higher density of available carbon to form 
this (simple) carbon-chain molecule in the higher mass-loss rate 
C stars. 
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Fig.B.l. New data from HIFI and IRAM for C stars not modelled in this paper, plotted with respect to LSR velocity. 







Fig.B.2. New data from HIFI and IRAM for S stars not modelled in this paper, plotted with respect to LSR velocity. 
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Fig.B.3. New data from HIFI and IRAM for M stars not modelled in this paper, plotted with respect to LSR velocity. 
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Fig. B.3. Cont. New data from HIFI and IRAM for M stars not modelled in this paper, plotted with respect to LSR velocity. 
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Fig.B.4. New SiO data from HIFI, plotted with respect to LSR velocity. 
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Table B.l. ObsIDs for Herschel observations. 


Source 

OhsID CO(5 ^ 4) 

OhsID CO(9 ^ 8) 

OhsID CO(14 ^ 13) 

VI968 Cyg 

1342210073 

1342221424 


AFGL 292 

1342200903 

1342213353 


V701 Cas 

1342201532 

1342227533 


V1259 Ori 

1342218903 

1342227530 


V688 Mon 

1342218902 

1342229929 


BK Vir 

1342210683 

1342212125 


BX Cam 

1342230363 

1342227532 


EP Aqr 

1342210079 

1342220501 


GXMon 

1342228562 

1342228602 


IRAS 06192+4657 

1342230362 

1342227531 


AI Vol 

1342204015 

1342200911 

1342210044 

1342229787 

II Lup 

1342202049 

1342227540 


VIIII Oph 

1342230369 

1342229924 


KU And 

1342200923 

1342213361 


V370 And 

1342201528 

1342213354 


V384 Per 

1342204006 

1342204007 

1342227534 

1342216333 


NV Aur 

1342217689 

1342218627 


GY Cam 

1342214328 

1342218626 


VI300 Aql 

1342210080 

1342216809 


L 2 Pup 

1342202055 

1342210049 

1342231773 

n Gru 

1342210084 

1342210043 

1342220504 

R And 

1342200924 

1342213362 


RCrt 

1342210087 

1342212121 


RCyg 

1342200920 

1342210034 


R Gem 

1342228561 

1342228601 


RHor 

1342200910 

1342213350 


RHya 

1342200913 

1342212120 

1342223432 

RLMi 

1342220512 

1342220498 


R Leo 

1342210684 

1342220493 


R Lep 

1342214326 

1342216330 


RScl 

1342210695 

1342200909 

1342213347 

1342221454 

RR Aql 

1342216806 

1342216808 


RS Cnc 

1342220513 

1342220497 


RT Vir 

1342200915 

1342212126 


RV Aqr 

1342210078 

1342218413 


RXBoo 

1342200916 

1342212129 

1342223433 

S Cas 

1342204009 

1342201533 

1342213358 


S CMi 

1342220514 

1342220496 


SCrB 

1342200917 

1342212128 


SW Vir 

1342200914 

1342212124 


TCep 

1342201535 

1342210666 


TXPsc 

1342210686 

1342222348 


UHya 

1342210088 

1342212122 


UU Aur 

1342230361 

1342229926 


V Aql 

1342230383 

1342229921 


VCrB 

1342200918 

1342214438 


V1943 Sgr 

1342216804 

1342216812 


V821 Her 

1342230368 

1342215900 


W And 

1342201527 

1342213355 


XHer 

1342200919 

1342210033 


XTrA 

1342202053 

1342216329 


YLyn 

1342230360 

1342229925 


YUma 

1342210067 

1342212127 
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Table B.3. HIFI and IRAM CN and H3CN line group observations 


Star 

I 1/2 —> O 1/2 
[Kkms-'l 

I 3/2 —> O 1/2 
[K km s“*] 

23/2 —> I 1/2 
[K km s“*] 

CN 

23/2 —> I 3/2 
[Kkm s“*] 

25/2 —> I 3/2 
[K km s“*] 

59/2 —» 47/2 

[K km s“'] 

5ii/2 —> 49/2 
[K km s“'] 

HC 3 N 

12 ^ 11 
[K km s“'] 

RScl 

15.6(1.0) 

13.2(1.0) 

31.4 (0.9) 

7.1 (0.9) 

25.2 (0.9) 

3.0 (0.1) 

1.7 (0.1) 

< 1.0 

V701 Cas 

3.7 (0.3) 

8.5 (0.3) 

5.0 (0.4) 

<0.4 

6.3 (0.4) 

< 0.2 

< 0.2 

1.76 (0.17) 

V384 Per 

21.6 (0.2) 

28.4 (0.2) 

21.9 (0.3) 

7.0 (0.3) 

23.3 (0.3) 

< 0.2 

< 0.2 

3.72 (0.12) 

GY Cam 

13.8 (0.2) 

10.9 (0.2) 

16.2 (0.3) 

3.7 (0.3) 

12.9 (0.3) 

< 0.2 

< 0.2 

0.69 (0.17) 

R Lep 

11.8 (0.4) 

11.7 (0.4) 

37.2 (0.6) 

6.8 (0.6) 

29.5 (0.6) 

2.0 (0.1) 

1.2 (0.1) 

<0.4 

VI259 Ori 

6.1 (0.2) 

12.6 (0.2) 

4.6 (0.3) 

<0.3 

8.9 (0.3) 

< 0.2 

< 0.2 

3.13(0.16) 

V688 Mon 

4.2 (0.3) 

8.7 (0.3) 

8.1 (0.2) 

<0.2 

11.6 (0.2) 

< 0.2 

< 0.2 

2.40 (0.15) 

UHya 

<0.4 

1.4 (0.4) 

17.1 (1.0) 

< 1.0 

11.6(1.0) 

1.8 (0.2) 

1.1 (0.2) 

<0.4 

XTrA 






1.0 (0.2) 

0.7 (0.2) 


V821 Her 

15.8 (0.5) 

22.6 (0.5) 

17.1 (0.7) 

<0.7 

25.2 (0.7) 

< 0.2 

< 0.2 

2.26 (0.24) 

V Aql 

4.4 (0.2) 

4.9 (0.2) 

10.6 (0.3) 

2.4 (0.3) 

13.0 (0.3) 

<0.1 

<0.1 

0.26 (0.08) 

VI968 Cyg 

3.0 (0.3) 

5.6 (0.3) 




< 0.2 

< 0.2 

1.06 (0.14) 

RV Aqr 

14.7 (0.2) 

15.7 (0.2) 

22.2 (0.4) 

4.7 (0.4) 

20.2 (0.4) 

<0.1 

<0.1 

0.28 (0.15) 

S Cas 

<0.3 

1.4 (0.3) 

11.8(0.5) 

2.4 (0.5) 

11.6 (0.5) 

< 0.2 

< 0.2 

<0.3 


Notes. The value in brackets after the flux gives the integrated noise RMS. An ellipsis (...) indicates that the line was not observed for the indicated 
star. An * indicates the line was detected but contaminated by an artefact and the integrated intensity cannot be relied upon. 
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Fig.B.5. New '^CO data from IRAM, plotted with respect to LSR velocity. 
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Fig.B. 6 . New SiS and HC 3 N data from IRAM, plotted with respect to LSR velocity. Note that the peak at- 80 km s * in the GY 

Cam H 3 CN spectrum is an artefact and not part of the H 3 CN line. 
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Fig. B. 7 . New CN data from HIFI and IRAM, plotted with respect to LSR velocity of the reddest component. In the case of the HIFI 
lines, both the (59/2 ^ 47/2) and ( 5 11/2 —> 49/2) lines are plotted together at the rest frequency of the former. 











































































































































































































































T. Danilovich et ah: Observations and models of circumstellar CO line emission of AGB stars in the SUCCESS programme, Online Materialp 18 









Velocity [km/s] 



Velocity [km/s] 


Velocity [km/s] 


Velocity [km/s] 


Velocity [km/s] 


Fig. B. 7 . Cont. New CN data from HIFI and IRAM, plotted with respect to LSR velocity of the reddest component. In the case of 
the HIFI lines, both the (59/2 ^ 47/2) and ( 5 ii /2 ^ 49/2) lines are plotted together at the rest frequency of the former. 
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Appendix C: Supplementary line data 

As discussed in Sect. | 2 . 3 | we included substantial archival data 
in our modelling procedure to find the models which best fit the 
widest range of data possible. The archival data we used to con¬ 
strain our models is listed in Table IC.ll 


Table B.2. IRAM and HIFl ^^CO, SiO and SiS line observations 


Star 

"CO (1 ^ 0) 
[K km s“'] 

SiO (13 ^ 12) 
[Kkm s“*] 

SiS (6 ^ 5) 
[K km s“‘] 

C Stars 

RScl 

4.7 (0.4) 

<0.1 

< 0.4 

V701 Cas 

* 2.7 (0.2) 

< 0.22 

0.94 (0.15) 

V384 Per 

3.9 (0.2) 

0.73 (0.18) 

0.77 (0.14) 

GY Cam 

2.0 (0.2) 

< 0.20 

0.38 (0.15) 

V1259 Ori 

3.5 (0.2) 

< 0.21 

1.40 (0.13) 

V688 Mon 

* 1.5 (0.2) 

<0.19 

0.41 (0.12) 

V821 Her 

7.7 (0.4) 

0.29 (0.16) 

1.33 (0.22) 

RV Aqr 

2.2 (0.2) 

<0.13 

0.39 (0.14) 

S Stars 

R And 

0.83 (0.29) 

<0.16 

<0.3 

S Cas 

1.3 (0.2) 

< 0.18 

< 0.2 

R Gem 

0.34 (0.26) 

<0.23 

<0.3 

RS Cnc 

1.0 (0.3) 

0.90 (0.16) 

<0.3 

RCyg 

0.41 (0.21) 

<0.16 

< 0.2 

M Stars 

KU And 

3.7 (0.2) 

<0.16 

0.55 (0.15) 

V370 And 


1.07 (0.22) 


AFGL 292 

0.16(0.10) 

<0.23 

<0.1 

RHor 


0.41 (0.18) 


NV Aur 

7.9 (0.2) 

<0.13 

1.25 (0.12) 

BX Cam 

2.9 (0.2) 

0.85 (0.22) 

0.40 (0.11) 

GXMon 

6.2 (0.2) 

< 0.18 

1.04(0.12) 

L 2 Pup 


1.13 (0.19) 


RLMi 


0.54 (0.22) 


R Leo 


1.46 (0.18) 


RCrt 


1.73 (0.17) 


BK Vir 

0.04(0.13) 

0.29 (0.17) 

<0.1 

YUMa 


0.10(0.18) 


RT Vir 


0.87 (0.18) 


SW Vir 


1.50 (0.22) 


RHya 


0.72 (0.18) 


RXBoo 


1.63 (0.20) 


XHer 


0.79 (0.27) 


Vllll Oph 

5.6 (0.3) 

<0.19 

<0.3 

V1943 Sgr 


0.37 (0.17) 


V1300 Aql 

11.0 (0.4) 

<0.15 

1.53 (0.22) 

TCep 


0.32 (0.17) 


EP Aqr 


2.00 (0.21) 



Notes. The value in brackets after the flux gives the integrated noise 
RMS. An ellipsis (...) indicates that the line was not observed for the 
indicated star. * indicates that flux has been corrected for ISM emission. 
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Table C.l. Archival data of other CO observations of the stars used in our modelling. 


Star Transition /^b Telescope 

[K km s“'] 

C Stars 


V384 Per 

1 ^0 

52.7 

IRAM 


1^0 

7.8* 

NRAO 


1 ^0 

25.5* 

OSO 


2 ^ 1 

83.6 

IRAM 


2 ^ 1 

61.9 

JCMT 


3^2 

32.8 

CSO 


3^2 

40.1 

JCMT 


4^3 

78.7 

JCMT 

GY Cam 

2 ^ 1 

14.3 

CSO 

R Lep 

1^0 

6.2 

SEST 


2^ 1 

27.6 

APEX 


2 ^ 1 

18.1 

SEST 


3^2 

32.3 

APEX 


4^3 

36.6 

APEX 

V1259 Ori 

2^ 1 

36.3 

APEX 


2 ^ 1 

16.5 

CSO 


3^2 

32.5 

APEX 


3^2 

32.2 

APEX 


4^3 

30.4 

APEX 

UU Aur 

1^0 

18.8 

IRAM 


1 ^0 

7.9 

OSO 


2 ^ 1 

15.0 

APEX 


2^ 1 

39.0 

IRAM 


3^2 

21.6 

APEX 


3^2 

22.7 

JCMT 

V688 Mon 

3^2 

15.4* 

CSO 

Al Vol 

1 ^0 

25.57 

SEST 


2^ 1 

56.5 

APEX 


2 ^ 1 

59.1 

APEX 


2 ^ 1 

50.00 

SEST 


3^2 

57.8 

APEX 


4^3 

59.1 

APEX 

UHya 

1 ^0 

5.4 

SEST 


2 ^ 1 

17.9 

APEX 


2^ 1 

48.8 

IRAM 


2^ 1 

20.2 

JCMT 


2 ^ 1 

13.8 

SEST 


3^2 

25.8 

APEX 


3^2 

27.2 

APEX 


3^2 

29.3 

APEX 


4^3 

23.1 

APEX 


4^3 

30.6 

APEX 


7^6 

25.1 

APEX 

XTrA 

1^0 

2.5 

SEST 


2 ^ 1 

10.9 

APEX 


2 ^ 1 

11.8 

SEST 


3^2 

15.0 

APEX 


3^2 

12.3 

APEX 


3^2 

15.3 

SEST 


4^3 

19.5 

APEX 

II Lup 

1^0 

61 

SEST 


2^ 1 

118.4 

APEX 


2 ^ 1 

128.12 

APEX 


2 ^ 1 

151 

SEST 


3^2 

163.0 

APEX 


3^2 

144.4 

APEX 


3^2 

144.4 

APEX 


3^2 

130.3 

APEX 


3^2 

122.71 

APEX 


3^2 

129.3 

SEST 


4^3 

155.9 

APEX 


4^3 

128.4 

APEX 


4^3 

147.0 

APEX 


6^5 

20.2 

HIFI 

* indicates uncorrected ISM contamination 


Reference 


Olofsson et al. i 1993 

1 

Schoier & Ulotsson 

2001 

Ulotsson et al. 

1993 

( 

Ulotsson et al. 

1993 

( 

Schdier & Ulot 

sson 

2001 

Knapp et al. i 1998 ( 

Schoier & Ulotsson 

2001 


hi. Olotsson, privat e communication 


Knapp et al. 

\D 

\D 

00 

Ulotsson et al. 

^1993 


Pointing Catalogue 


19931 


Olofsson et al. 
Pointing Catalogue 
Pointing Catalogue 
Pointing Catalogue 


Knapp et al. 1 1998 1 
Pointing Catalogue 
Archive 


Olofsson et al. 

1993 

Ulotsson et al. 

1993 


Pointing Catalogue 


Olofsson et al. 


119931 


Pointing Catalogue 


Knapp et al. 

1998 

Woods et al. 

2003 


Hans pointing data 


Woods et al. 1 2003 1 
Pointing Catalogue 
Pointing Catalogue 


Olofsson et al. 


19931 


Pointing Catalogue 


Schdier & Olofsson 

2001 

Schdier & Olotsson 

2001 

Schdier & Ulotsson 

2001 


Archive 


De Beck et al. 


(20101 


Pointing Catalogi 


De Beck et al. 

20101 

De Beck et al. 

20101 

Ulotsson et al. 

1993 

Pointing Catalc 

)gue 


19931 


Olofsson et al. 
Pointing Catalogue 
Archive 


Schdier & Olofsson 
Pointing Catalogu e 


20011 


Ryde et al. 1 1999 1 
Pointing Catalogue 
Archive 


Ryde et al. i 1999 

( 

Ramstedt & Ulotssor 

2014 

De Beck et al. 

2I7TU 


De Beck et al. 

2010 


Archive 


Ryde et al. 

1999 

( 

De Beck et al. 

2010 

De Beck et al. 

2010 


Herschel Science Archive (HIFISTARS) 
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Table C.l. continued... 


Star 


Transition 


fmb 
[K km s“'] 


Telescope Reference 


VCrB 


V821 Her 


V Aql 


VI968 Cyg 
RV Aqr 


5 Stars 
R And 


S Cas 

W And 
R Gem 


YLyn 
RCyg 
M Stars 


16- 
1 ■ 
1 ■ 
2 ■ 
2 - 
3 ■ 
1 ■ 
2 ■ 
2 - 
3 ■ 
3 ■ 
3 ■ 
4- 
4- 
4- 
7 ■ 
1 ■ 
1 ■ 
2 - 
2 ■ 
2 ■ 
3 ■ 
3 ■ 
1 ■ 
2 ■ 
1 ■ 
2 - 
3 ■ 

2 - 
3 ■ 
4- 
1 ■ 
3 ■ 
1 ■ 
1 ■ 
2 ■ 
3 ■ 
1 ■ 
1 ■ 


• 15 
•0 
•0 

• 1 
• 1 
• 2 
•0 
• 1 
• 1 
• 2 
• 2 
• 2 

• 3 

• 3 

• 3 

• 6 
•0 
•0 
• 1 
• 1 
• 1 
• 2 
• 2 
•0 
• 1 
•0 
• 1 
• 2 

• 1 
• 2 

• 3 
•0 

• 2 
•0 
•0 
• 1 
• 2 
•0 
•0 


7^6 

151.43 

APEX 

De Beck et al. 

K) 

O 

O 

1 

7^6 

82.1 

APEX 

De Beck et at. 

21I[U 

( 


17.2 

11.7 

2.7 

5.4 

18.4 

10.9 

84.1 

52.1 

93.4 

60.3 

56.7 

55.9 

66.1 
66.0 

81.3 

49.8 

3.2 

2.8 
8.8 
9.0 
8.1 

11.5 
9 

46.4 

88.6 
7.5 
18.1 
18.6 

32.0 

43.0 

25.6 

13.8 
31.0 

3.8 

2.4 

4.5 

5.2 
4.1 
4.4 


HIFI 

IRAM 

OSO 

APEX 

IRAM 

APEX 

IRAM 

APEX 

IRAM 

APEX 

APEX 

APEX 

APEX 

APEX 

APEX 

APEX 

OSO 

SEST 

APEX 

JCMT 

SEST 

APEX 

JCMT 

IRAM 

IRAM 

SEST 

SEST 

SEST 

JCMT 

JCMT 

APEX 

OSO 

JCMT 

OSO 

OSO 

APEX 

APEX 

OSO 

OSO 


Herschel Science Archive (HIFISTARS) 


Schdier & Olofsson 

2001 

Ototsson et al. 

1993 

( 


Olofsson et al. 
Archive 


119931 


Neri et al. 


(1998 


Pointing Catalog i 


Neri et al. 

1998 

1 

De Beck e 

t al. 

i2 

010 


Archive 


(20101 


De Beck et al. 
Pointing Catalogue 
Archive 


De Beck et al. 

20101 

Olofsson et al. 

1993 

Ototsson et al. 

1993 


Pointing Catalogue 


Schdier & Olofsson 

2001 

Olofsson et al. 

1993 

( 


Pointing Catalogue 


Schdier & Olofsson 

2001 

Neri et al. 

TW5 

3 


Neri et al. 

199i 

3 

Ototsson et al. 

1993 

1 

Olotsson et al. 

1993 

1 

Schdier & Olot 

5 son 

2001 


Ramstedt et al. 

2009 

Ramstedt et al. 

2009 


Ramstedt et al. 

2009 

Ramstedt et al. 

2009 

Ramstedt et al. 

2009 

Ramstedt et al. 

2009 


Pointing Catalogue 


Ramstedt et al. 

2009 

Ramstedt et al. 

2009 


KU And 

1^0 

24.8 

OSO 


2^ 1 

14.0 

CSO 


2 ^ 1 

41.5 

JCMT 

RHor 

2 ^ 1 

19.0 

APEX 


3^2 

30.8 

APEX 


3^2 

23 

CSO 


4^3 

33.5 

APEX 

NV Aur 

2 ^ 1 

37 

JCMT 


2^ 1 

36.8 

JCMT 


2^ 1 

36.3 

JCMT 


3^2 

39 

JCMT 


3^2 

35.7 

JCMT 


3^2 

38.8 

JCMT 


4^3 

35.4 

JCMT 

BX Cam 

1 ^0 

20.0 

OSO 


2 ^ 1 

25.5 

JCMT 


3^2 

39.4 

JCMT 

GXMon 

1 ^0 

31 

OSO 


2^ 1 

56.7 

APEX 


2 ^ 1 

52.9 

APEX 


2 ^ 1 

22.9 

CSO 


* indicates uncorrected ISM contamination 


Gonzalez Delgado et al. 

2003 

Knapp et al. 


Gonzalez De 

gado et al. 

2003 


Pointing Cata logue 


Young(1995i 


Kemper et al. i 

20031 


Gonzalez Delg 

ado et al. 

2003 

De Beck et al. 

217TU 



Kemper et al. 

2003 1 

Gonzalez Delg 

ado et al. 

2003 

De Beck et al. 

2010 



De Beck et al. 

2010 


Ototsson et al. 

1998 



Ramstedt & O 

otsson 

2014 

1 

Ramstedt & Ototsson 

2014 

Ramstedt et al 

ITOOK 




Archive 


Knapp et al. (19981 
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Table C.l. continued... 


Star 

Transition 

^mb 

[K km s“^] 

Telescope 


2^ 1 

61 

JCMT 


3^2 

64.5 

APEX 


3^2 

55.0 

APEX 


3^2 

71 

JCMT 


4^3 

54.3 

APEX 


4^3 

149 

JCMT 

Lt Pup 

1 ^0 

0.24 

SEST 


2^ 1 

3.9 

APEX 


2^ 1 

3.1 

JCMT 


2 ^ 1 

3.7 

SEST 


3^2 

00 

bo 

APEX 


3^2 

16.6 

JCMT 


4^3 

11.5 

APEX 

RLMi 

2 ^ 1 

6.1 

APEX 


2 ^ 1 

2.72 

CSO 


3^2 

9.1 

APEX 


3^2 

12.9 

APEX 

R Leo 

1^0 

4.1 

IRAM 


1 ^0 

2.1 

OSO 


2 ^ 1 

14.8 

APEX 


2^ 1 

28.5 

IRAM 


2^ 1 

15.0 

JCMT 


3^2 

30.3 

APEX 


3^2 

32.7 

APEX 


3^2 

21.9 

CSO 


3^2 

41.6 

JCMT 


4^3 

39.5 

APEX 


4^3 

28.1 

CSO 


6^5 

31.0 

CSO 

SCrB 

2^ 1 

6.4 

APEX 


2 ^ 1 

2.53 

CSO 


2 ^ 1 

12.7 

IRAM 


3^2 

10.6 

APEX 


3^2 

10.5 

APEX 


3^2 

9.1 

CSO 


4^3 

4.9 

CSO 

Vllll Oph 

1^0 

21.5 

OSO 


2^ 1 

30.6 

APEX 


2 ^ 1 

42.0 

JCMT 


3^2 

40.5 

APEX 


3^2 

59.1 

JCMT 


4^3 

46.1 

APEX 

RR Aql 

2 ^ 1 

12.5 

APEX 


2 ^ 1 

16.6 

JCMT 


3^2 

16.1 

APEX 


3^2 

23.5 

JCMT 


4^3 

14.3 

APEX 

V1943 Sgr 

1 ^0 

1.25 

SEST 


2^ 1 

5.17 

SEST 

VI300 Aql 

1^0 

16 

OSO 


2 ^ 1 

38.1 

APEX 


2 ^ 1 

45 

JCMT 


3^2 

37.0 

APEX 


3^2 

35.3 

APEX 


3^2 

35.9 

APEX 


3^2 

27 

JCMT 


4^3 

33.5 

APEX 


4^3 

34.8 

APEX 


4^3 

34.5 

APEX 


4^3 

22 

JCMT 


7^6 

19.7 

APEX 

TCep 

1^0 

1.7 

OSO 


* indicates uncorrected ISM contamination 


Reference 


(2008 


Ramstedt et al. 
Pointing Catalogue 
Archive 


Ramstedt et al. 


(2008 


Pointing Catalogue 


Ramstedt et al. 

2008 

1 

Kerschbaum et at. 


96 


H. Olofsson, private c ommunication 


Olofsson et al. 
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20021 


120021 


Olofsson et al. 
Pointing Catalogue 
Pointing Catalogue 


Knapp et al. 1 1998 1 
Pointing Catalogue 
Archive 


Tevssier et al. 

200611 

Gonzalez Delgado et al. 

bJ 

o 

o 

OJ 


Gonzalez De 

Igado 

et al. 

o 

o 

Pointing Catalogue 
Archive 


Knapp et al. 
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Gonzalez De 

gado et al. 

2003 


Pointing Catalogu'e 


Young 1 1995 


leyssier et al. 

(2006 

Pointing Cata 

loeue 

Knapp et al. 1 

19981 

Nenetal.lTggS 


Archive 


Young 

1995 

1 

Young 

1995 

1 


Gonza 

ez De 

gado et al. 

2003 


Pointing Catalogu'e 


Gonzalez Delgado et al. (|2003 i 
Pointing Catalogu'e 


20141 


Ramstedt & Olofsson 
Pointing Catalogue 
Pointing Catalogue 
H. Olofsson, private communication 
Pointing Catalogue 
H. Olofsson, private communication 


Kerschbaum & Olofsson 

1999 

Kerschbaum & 

Olotsson 

1999 

Ramstedt et al. 



Ramstedt et al. 

(2008 

De Bectc et al. 

2010 1 


Archive 


Ramstedt et al. 

(2008 

De Bectc et al. 

2010 1 

Pointing Catalogue 
Archive 

Ramstedt et al. 

(20081 

De Bectc et al. 

20Tnr 

Olotsson et al. 

(1998# 




























































































































































